Strength distributions and symmetry breaking in coupled 

microwave billiards 

B. Dietz/ T. Guhr,^ H.L. Harney,^ and A. Richter^'Q 

^Institut fiir Kernphysik, Technische Universitdt Darmstadt, D-64289 Darmstadt, Germany 
^ Matematisk Fysik, LTH, Lunds Universitet, S-22100 Lund, Sweden 
^ Max-Planck-Institut fiir Kernphysik, D- 69029 Heidelberg, Germany 

(Dated: February 8, 2008) 

Abstract 

Flat microwave cavities can be used to experimentally simulate quantum mechanical systems. 
By coupling two such cavities, we study the equivalent to the symmetry breaking in quantum 
mechanics. As the coupling is tunable, we can measure resonance strength distributions as a 
function of the symmetry breaking. We analyze the data employing a qualitative model based on 
Random Matrix Theory (RMT) and show that the results derived from the strength distribution 
are consistent with those previously obtained from spectral statistics. 
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The breaking of quantum mechanical symmetries represents a prominent object of study 
in a rich variety of systems, ranging from high energy to condensed matter physics. Some- 



times it is possible to determine the size of the symmetry breaking by analyzing statistical 
fvables. We mention parityviolation ll , t 

m 

and isospin mixing |5|, |6|, |7|, |8|, |9| 



observables. We mention parityviolation ll| , the breaking of atomic and molecular symme- 
tries and isospin mixing jj, O, 0, Q, mm nucle . Symmetry breaking influences the 
spectral statistics as well as the wave function statistics, such as the distributions of partial 
widths and transition matrix elements. While symmetry breaking cannot be controlled or 
tuned in nuclei, atoms and molecules, this was possible for elastomechanical resonances in 
quartz crystals by successively breaking crystal symmetries An analo gy t o the nuclear 



case of symmetry breaking was given by two coupled microwave cavities [ll| with a tun- 
able coupling. The statistical properties of both systems are fully equivalent to those of a 
quantum system [3, ll^ • Wheras these two investigations focused on the spectral statistics, 
width distributions for two different types of resonances in elastic aluminum plates were 
studied in Ref. for a varying degree of mixing. 

Here, we present experimental results on the distribution of products of partial widths 
in two coupled chaotic microwave cavities. First, we measure the distribution for different 
couplings. It can be normalized such that it only depends on the symmetry breaking. As 
this is different from Ref. jl^, we can apply a qualitative statistical model which extends 



the model of Ref. [8|] in order to, second, extract the size of the symmetry breaking from 
the data. We then show, third, that the symmetry breaking thereby obtained is consistent 
with the one previously found from the spectral statistics. 

In the experiment, we used two flat cylindrical microwave cavities having the shape of 
the quarter of a Bunimovich stadium ^llj. In both resonators (see Fig. ^ the radius of 
the quarter circle is 0.2 m. The ratios of the length of the rectangular part to the radius 
are 71 = 1 and 72 = 1.8, respectively. Therefore the level densities increase with different 
slopes as functions of frequency [3, ll^- The cavities were put on top of each other and 
circular holes, 4 mm in diameter, were drilled through the walls of both resonators. The 
coupling was realized by a niobium pin, 2 mm in diameter, penetrating through the holes 
into both resonators. Due to the ring-shaped gap between the niobium pin and the hole 
surface, the pin acts like an antenna, which in the experimental frequency range supports 
one TEM-mode. The coupling is controlled by the penetration depth 

Information about the partial widths is obtained from transmission spectra. At a given 
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frequency /, the relative power transmitted from antenna a to antenna h is proportional to 
the absolute square of the scattering matrix element, Pout,b/-Pin,a ~ |'S'afe(/)p. For sufficiently 
isolated resonances, one has 



for / close to the frequency of the /i-th resonance. The quantities and F^^ are 
;he partial widths related to the antennae a and 6, F^ is the total width of the resonance 
ItI i. In the experiment, three antennae were attached to each resonator, where half of the 
microwave power was fed into each resonator (see insert of Fig. E}. Thereby, altogether 
six transmission spectra were obtained for 4 different couplings; in the order of increasing 



coupling these are denoted as (8,0), (5,3), (4,4) and (5,8) in Ref. Here, the pair (xi,X2) 
denotes the penetration depth in mm of the pin into either cavity. The transmission spectra 
were measured up to a maximum frequency of 17.5 GHz. The resonance strengths F^^F^fe, 



i.e. the products of the partial widths, are determined as described in 18|. Resonances with 
peak heights below a certain value may not be detected, implying that some strength is 
missing in the tails of the distributions We will show that this has no effect on the 

results. 

We work with the resonance strengths F^^F^fe and their distributions. They yield the 
same information as the distributions of the partial widths F^^ themselves. The data are 
unfolded as in Refs. 
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To interpret the distdbufon of the emp.dea. data F^,. we employ a statistiea. model 

which is a special case of the Rosenzweig-Porter model [2]. A symmetry is associated with 
a quantum number. If it suffices to consider only two different values of it, the Hamiltonian 
can be written in the form 

(Hi Q \ / r\ 

where the first part preserves the symmetry and the second part breaks it with the parameter 
a. For isospin mixing in nuclei. Hi and H2 correspond to the sub-Hamiltonians for two 
isospin quantum numbers, while V accounts for the Coulomb interaction that mixes isospin. 
In our case the "symmetry" preserving situation is simply realized by considering the states 
of each uncoupled cavity as eigenstates of a "symmetry operator" . Thus, Hi and H2 model 
the dynamics in the two cavities, without the coupling. We choose their dimensions A^^i and 



N2 different because the level densities in the cavities have different slopes as a function 
of frequency. For later purposes, we introduce the fractional densities gi = Ni/{Ni + A^2) 
and g2 = A^2/(^i + ^2)- As the whole system is time-reversal invariant and the dynamics 
in the Bunimovich billiards is fully chaotic, we represent Hi and H2 by real -sy mmetric 



random matrices drawn from the Gaussian Orthogonal Ensemble (GOE) [12, ll9|, |20|, |21 |. 
The coupling is modeled by the off-diagonal blocks V and V'^ in Eq. (0), where the matrix 

V is real with no symmetries and has dimension Ni x N2. In the experiment the modes in 
one resonator are coupled to those in the other via one TEM-mode. We model the coupling 
of the A^^i and N2 resonator modes with this TEM mode by an A'^i-dimensional vector v 
and an A^2- dimensional vector w. The coupling matrix V then acquires the dyadic structure 

V = vw'^. The matrix V has rank M = 1. The entries of v and w are chosen as Gaussian 
random numbers. The distribution of the elements Vnm is 

q{VnmW) dVnm = ^ dVnm , (3) 



where Kq is the modified Bessel function of zeroth order |22]. The variance a of the elements 
of V is adjusted exactly as described in jsl, O]. The vertical bar on the l.h.s. separates the 
argument of q into the random variable Vnm and the parameters, i.e. in this case a. Universal 
spectral fluctuations have to be measured on the scale of the local mean level spacing D. 
Accordin gly, the parameter measuring the size of the coupling or the symmetry breakingis 
\ = a/D l23|. It connects to the spreading width via the relation r-*- = 27ra'^/D = 2'k\'^D l5|. 

nn 

Spectral fluctuations have been calculated perturbatively in A in Refs. |24L I25l | and exactly 
for the time-reversal non-invariant case in Refs. 0, |2q. For the analysis of the spectral 



properties of the coupled microwave billiards 



ll| , a model Hamiltonian of the form Eq. (j21) 



with a coupling matrix V of rank A^^i was applied. While, in the relevant range of A-values 
the spectral properties do not depend on the rank of V , the statistical properties of the 
eigenvector components deviate for M <^ Ni. ^ 

Accordingly we apply and extend the qualitative model of Ref. |S] where the statistics of 
transition strengths was considered. In a fully chaotic system, the partial width distribution 
converges to the Porter-Thomas form |l2 




PT(tJr,)dt„ = , exp (4) 



for a large level number. We write ta for the partial width instead of F^a, because the 
distribution does not depend on the resonance index \i. In practice the distribution is 



obtained from the sample of all resonances. Its first moment equals r^. For two coupled 
chaotic systems the partial width distribution p{ta\X,Ta) involves the symmetry breaking 
Hamiltonian in Eq. Q for large level numbers. It is easily constructed in the limiting case 
without symmetry breaking; i.e. for A = we have 



The distribution in Eq. (0) again has the expectation value r^. We now model the case of 
small symmetry breaking A by an interpolating ansatz: We expect the Porter-Thomas dis- 
tributions in Eq. (jSj) to maintain their shape, such that only their width parameters change. 
The delta function in Eq. (0) acquires a width and develops into a non-singular function. Nu- 
merical studies have lead us to approximate it by ^8J Po(^|p) dt = Kq {\Jt/ p^/ i'n'y/tp) dt, 
where p stands for the variance parameter. We arrive at the model 



The shape of the whole distribution is determined by the quantities ni{f]), ^2(^7) and rj = 
77(A) which all depend on A. With the help of the limiting case Eq. (0), we can come 
up with educated guesses for these functions. We must have 77(0) = and furthermore 
Kj(0) = I/qj, j = 1,2. As the functions Kj ought to be even in rj, we choose Kj{ri) = 
l/gj + {1 — l/gj) rj"^ , j = 1, 2. To construct the function rj = //(A), we fit the ansatz given 
in Eq. (jH)) including these choices to Monte Carlo simulations of the distributions involving 
the Hamiltonian in Eq. Q for A^^i + A'"2 = 300. We do that for 21 values of A in the interval 
[0.05, 0.25] which is the relevant parameter range in the experiment. The fractional densities 
in the experiment are gi = 0.59 and g2 = 0.41. In Fig. |21we display the values thus obtained 
for //(A). The functional form is well described by the polynomial 1]{X) = 2.57A — 1.98A^. 

From the distributions for two partial widths ta and tb, say, we obtain the distribution for 
the resonance strength y = tatb as in by a filter integration. As we know the fractional 
densities gi and g2 in the experiment, we arrive at a qualitative model for the resonance 
strength distribution that depends only on one single parameter — the measure A for the 
symmetry breaking. The distributions in Eq. © and accordingly the strength distribution 




(5) 



+2gig2Po(ta\Tar]^) . 



(6) 
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diverge for — > 0. Thus, for their graphical representation we use the logarithmic variable 

Despite its simplicity, the qualitative model in Eq. (jH)) yields a satisfactory description. 
We demonstrate this in Fig. |21 where a Monte Carlo simulation of resonance strength distri- 
butions obtained from the random matrix Hamiltonian in Eq. Q with the same fractional 
densities as in the experiment is compared to the calculated strength distribution for four 
different values of A. The strongest deviations between both curves are observed for the 
smallest value of A. For values A < 0.03 the description by the qualitative model ceases to 
be satisfactory. But, most importantly, the position of the maximum is described very well 
in all four examples. It is particularly this feature which makes the qualitative model in 
Eq. © useful for the analysis of experimental data. The fit of the model to the experimental 
data is displayed in Fig. EJ Due to the above mentioned missing strength in the left tails, 
the shape of the distribution is not described quantitatively. In fact, for the determination 
of A only strengths in a 2; interval [-3,1.5] were taken into account, where the probability of 
missing strength is small. The corresponding experimental distributions are shown together 
with the RMT model fits in the insets of Fig. E]b)-d). For very small values of A (cf. Fig. |3] 
a)) the distributions agree fairly well in this interval. We emphasize again that the peak 
positions in the chosen interval of z values uniquely determine the coupling strengths. We 
thus expect reliable results for the parameter A. To carefully check how much the missing 
strength influences the extracted values of A, we amended the qualitative model by taking 
care of the experimental thresholds using an analytic ansatz 3| • As expected, the influence 
of the thresholds on the extracted values of A turned out to be negligible, except when A 
is essentially zero. As only the position of the maximum in the strength distribution is 
relevant, we have some freedom in choosing the distribution Pq in Eq- ©• We also tried a 
Porter-Thomas distribution, but the KQ-one describes the shape of the resonance strength 
distribution better. 



We finally compare the coupling strengths A 



tracted from the spectral correlations in Ref. 
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bund in the present work with those ex- 
22|. The values are given in Tab. HI The 



A values obtained in the present work are averaged over all six antennae combinations. The 
two analyses agree within the experimental errors. For the weakest coupling only an upper 
limit can be given. This is due to the fact that for essentially zero coupling a consider- 
able share of strength lies below the experimental threshold of detection. The consistency 
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for the other couphngs is an encouraging corrobaration of the analysis carried out in this 
contribution js^. 

In conclusion, we have measured resonance strength distributions for two coupled mi- 
crowave billiards modeling quantum systems with symmetry breaking. We analyzed the 
data with a qualitative model which depends only on one single free parameter — the size 
of the symmetry breaking. Our results are interesting from an additional point of view. The 
spectral correlations are more strongly affected by missing levels than partial widths and 
transition or resonance strength distributions are affected by missing strength. This is so 
because the latter do not comprise correlations, they are just densities. Thus the observ- 
ables related to the wave functions may provide more reliable information than the spectral 
correlations. We have shown, with data much better than are usually available, that the 
empirical information extracted from the wave function observables is consistent with that 
obtained from the spectral correlations. 
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FIG. 1: The two cavities are coupled with a tunable coupling 
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FIG. 2: Size of the tunable coupling A between the two cavities sketched in Fig. ^ The function 
r/(A) enters Eq. 
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FIG. 3: Monte Carlo simulation of resonance strength distributions (solid lines) compared to 
the calculated strength distribution (dashed lines) for the couplings A = 0.03 (a), A = 0.10 (b), 
A = 0.17 (c), A = 0.24 (d). 
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FIG. 4: The experimental resonance strength distributions (histograms) for the couphngs (8,0) 
(a), (5,3) (b), (4,4) (c) and (5,8) (d) fitted to the calculated strength distribution. The symmetry 
breaking parameter A equals A = in (a), A = 0.110 in (b), A = 0.125 in (c), and A = 0.185 in (d). 
The insets in b), c), and d) show experimental strength distributions together with RMT model 
fits for a z interval [-3,1.5] in which the probability of missing strength is small. 
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TABLE I: Symmetry breaking parameter A from the spectral correlations and from the resonance 



strength distributions. 



Physical 


A 


A 


coupling 


Ref. [nj 


Present work 


(8,0) 


< 0.029 


< 0.003 


(5,3) 


0.105 ±0.008 


0.116 ±0.003 


(4,4) 


0.130 ±0.007 


0.122 ±0.004 


(5,8) 


0.180 ±0.006 


0.195 ±0.007 
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